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Abstract. We compute the mass loss rate from a disk orbiting a stellar mass 
black hole assuming the flow is guided along magnetic field lines attached to the 
disk. We then estimate the Lorentz factor V of the wind at infinity. We find that 
r can reach high values only if severe constraints on the field geometry and the 
conditions of energy injection are satisfied. We discuss our results in the context of 
different scenarios for gamma-ray bursts. We mention a risk of "contamination" of 
the Blandford-Znajek process by wind material emitted from the disk. 



1. Introduction 

Among the sources which have been proposed to explain cosmic gamma-ray bursts 
(GRBs) the most popular are mergers of compact objects (neutron star binaries or 
neutron star - black hole systems) or massive stars which collapse to a black hole 
(hypernovae) jj], ||, ||]. In all cases, the resulting configuration is a stellar mass black 
hole surrounded by a thick torus made of stellar debris or of infalling stellar material 
partially supported by centrifugal forces. 

If black hole + thick disk configurations are indeed at the origin of GRBs the 
released energy will ultimately come from the accretion of disk material by the black 
hole or from the rotational energy of the hole itself extracted by the Blandford-Znajek 
mechanism. In a first step the energy must be injected into a relativistic wind. The 
second step consists in the conversion of a fraction of the wind kinetic energy into 
gamma-rays via the formation of shocks, probably inside the wind itself Q |J. In the 
last step the wind is decelerated when it interacts with the interstellar medium and 
the resulting (external) shock is responsible for the afterglow observed in the X-ray, 
optical and radio bands || . 

The origin of the relativistic wind is certainly the more complex of the three steps. 
A few possible ideas have been proposed but none is presently fully conclusive. If 
the burst energy comes from matter accretion by the black hole, the annihilation 
of neutrino-antineutrino pairs emitted by the hot disk could be a way to inject 
energy along the system axis, in a region which can be expected to be essentially 
baryon free due to the effect of centrifugal forces. The low efficiency of this process 
however requires high neutrino luminosities and therefore short accretion time scales 
. Another possibility is to suppose that disk energy is extracted by a magnetic field 
amplified by differential rotation to very large values (B ~ 10 15 G). A magnetically 
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driven wind could then be emitted from the disk with a fraction of the Poynting flux 
being eventually transferred to matter. An alternative to accretion energy could be 
to directly tap into the rotational energy of the black hole via the Blandford-Znajek 
mechanism. The available power then depends on the rotation parameter a of the black 
hole and on the intensity of the magnetic field pervading the horizon || . The purpose 
of this paper is to present an exploratory study of the case where a magnetically driven 
wind is emitted by the disk. Our approach will be extremely simplified in comparison 
to the complexity of the real problem so that our conclusions will have to be considered 
as indicative only. We nevertheless expect that we can identify the key parameters 
which control the baryonic load of such a wind and put constraints on the final values 
of the Lorentz factor which can be obtained. 



2. Dynamics of the wind from the disk to the sonic point 

To compute the mass loss rate and therefore estimate the amount of baryonic pollution 
we only need to follow the wind dynamics from the disk up to the sonic point. We 
write the wind equations with a number of simplifying assumptions: i) we assume 
that the disk is thin and that the field is poloidal with the most simple geometry, i.e. 
straight lines making an angle 9(r) with the plane of the disk, r being the distance 
from the foot of the line to the disk axis (Fig. 1). The flow of matter is then guided 
along the magnetic field lines; ii) we use non relativistic equations since even at the 
sonic point v s /c < 0.1 but we adopt the Paczynski-Wiita potential for the black hole; 
Hi) we consider that a stationary regime has been reached in the wind. 




Figure 1. Disk and field geometry 

We then write the three flow equations in a frame corotating with the foot of the 



line: 

• Conservation of mass 

pvs(y) = m 

• Euler equation 
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where y = Ijr and i is the distance along the field line; e is the specific internal energy, 
j(y) the total acceleration (gravitational + centrifugal) and q(y) the power deposited 
in the wind per unit mass. Different sources of heating can be present such as neutrino 
captures on nucleons (if the disk is optically thick to neutrinos) , neutrino-antineutrino 
annihilation and dissipation of kinetic or magnetic energy. Because the field and 
stream lines are coincident the function s(y) is easily related to the field geometry 
through the conservation of magnetic flux. Finally, m is the mass loss rate from 
the disk per unit surface. Our equation of state which includes nucleons, relativistic 
electrons and positrons and photons is computed from the expressions given by |^] . 

As long as the inclination angle 9 remains larger than 0\ ~ 60° (#1 is exactly 
60° if a newtonian instead of a Paczynski-Wiita potential is used for the black hole) 
the acceleration j(y) is negative up to y = y\ after which the centrifugal force 
dominates. The sonic point of the flow is located just below y\ (the relative difference 
(yi — y s )/yi never exceeds 1%). To solve the wind equations we first fix trial values of 
the temperature and density T s and p s at the sonic point from which we get v = v s - 
The position of the sonic point is obtained from the condition that the solution remains 
regular at y = y s . The mass loss rate m is then fixed and the inward integration along 
a field line can be started. We observe that at some position y = y cl -n the velocity 
begins to fall off rapidly while the temperature reaches a maximum T max < Td(^), 
where Tj}(r) is the disk temperature. We adjust the values of T s and p s with the 
requirement that y CI1 t should be as close as possible to and T max to Tn(r). 

The results presented below assume that the disk is optically thick to neutrinos 
between i?j n = 3 r g and i? m ax = 10 r g . No term for kinetic or magnetic energy 
dissipation have been included so that q(y) is limited to neutrino processes: capture on 
free nucleons, scattering on electrons and positrons, neutrino-antineutrino annihilation 
(heating) and neutrino emission by nucleons, annihilation of electron-positron pairs 
(cooling). The assumption that the disk is optically thick to neutrinos is probably 
justified for NS + NS or NS + BH mergers. It is much more questionable in the 
hypernova scenario, except for very high accretion rate or low values of the viscosity 
parameter (a < 0.01) as shown from the disk models computed by p0| . Our detailed 
results therefore only concern a specific case but we also obtain below a simple and 
general rough estimate of the mass loss rate for any kind of heating mechanism. 

The adopted temperature distribution Tb(r) corresponds to a geometrically thin, 
optically thick disk 



where T» is the temperature at r = r*. The mass of the black hole is Mbh = 2.5 M Q . 
3. The mass loss rate 

The solution for mass loss rate as a function of r, To(r) and 9(r) takes the form 
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where m 13 is the mass loss rate in units of 10 13 g.cm _2 .s _1 , x — r/r g and /ibh = 
A/bh/2.5 Mq. The geometrical function f[x, 0(x)] is normalized in such a way that it 
is equal to unity for x = 4 and 8 = 85°. The mass loss rate is extremely sensitive to 
the value of the disk temperature. The tenth power dependence is in agreement with 
what is found for neutrino driven winds in spherical geometry The dependence 
of to on inclination angle is also very strong as shown in Fig. 2 where to is represented 
(with T„ = 2 MeV and r* = A r g ) for two geometries of the field lines: constant 
6 = 85° and 9 decreasing from 90° to 80° between r — 3 r g and r = 10 r g . 
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Figure 2. Mass loss rate from the disk for a constant and decreasing inclination 
angle of the field lines. The disk temperature is 2 MeV at r = 4 r g 



Since additional sources of heating can be present in the wind (viscous dissipation, 
reconnection of field lines, etc) we have also obtained a very simple and general 
analytical expression for the mass loss rate j^j 

rh « -^g (6) 

where e is the rate of thermal energy deposition (in erg.cm _2 .s _1 ) between the plane 
of the disk (y = 0) and the sonic point at y s ~ yy\ A$ = $i — $o is the difference of 
potential (gravitational + centrifugal) between y = and y = y\. The g factor, which 
is of the order of unity, depends on the distribution of energy injection between y = 
and y = y 1 . 



4. Average Lorentz factor of the wind 

To estimate the Lorentz factor which can be reached by the wind one must be able 
to relate the injected energy to the mass loss rate. This can be done in the following 



■5 



way: we suppose that we observe a burst power in gamma-rays 
10 51 

47T 

Then, the power injected into the wind was 



£y = -73- £51 erg.s x .sr 1 (7) 



fO.l 
£0.06 



£ = 2 10 51 ^e 51 erg.s" 1 (8) 



where /q 1 and /° 05 are respectively the fraction of solid angle covered by the wind 
(in unit of 0.1) and the efficiency for the conversion of kinetic energy into gamma-rays 
(in unit of 0.05). Accretion by the black hole powers the wind but at the same time 
viscous dissipation heats the disk which cools by the emission of neutrinos. If neutrino 
losses represent a fraction a of the energy E injected into the wind we have (for an 
optically thick disk) 

fO.l /Tout n 

E v =aE = 2 10 51 j^ ae 51 = 2 -aT* (r) 2nrdr (9) 



■>1 

Substituting in (9) Eq. (4) for Tu (r) we obtain for the temperature at = 4 r 
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' = 1 72 Mbh /2 ( Ms<*<*i) ' MeV (10) 

v J 7 



The value of being known m can be computed as a function of of r for a given field 
geometry. The total mass loss rate from the disk is then 

/rout 
m(r) 2irrdr 

(ID 



2 -6 10 2 VIh ( ) ^ 



T. V 



2 MeV y 



where 



fr out /r g 

= / /[o;,0(a;)]a;da; (12) 

J nn/rg 

is a function of the field geometry. The average Lorentz factor is finally given by 

f-A-^V e-^a-^( f -^-\ m (13) 

T - mc*- t Mbh£51 a l^/r J (13) 

The value of T is 56 for a constant inclination angle = 85° and 250 if decreases 
from 90° to 80° between r — 3 and 10 r g . These numbers shows that large Lorentz 
factors can be reached but only under quite restrictive conditions: quasi vertical field 
lines, low a values, i.e. good efficiency for energy injection into the wind with little 
dissipation and necessity of beaming. 

More generally, Eq. (6) can also provide a simple and useful constraint on the 
terminal Lorentz factor. If the power e deposited below the sonic point represents a 
fraction x of the total power etot which is finally injected into the wind we get 

L«^|«^ (14) 



() 



Considering a line anchored at r = 4 r g with an inclination angle 9 = 85° we obtain 
y% = 2.182 and A$/c 2 = 0.18 which implies that x should not exceed 10 -3 to have 
r > 100! This is clearly a very strong constraint on any mechanism of energy injection. 
The wind however remains relativistic for x < 0.1 but its Lorentz factor is then much 
too low to produce a cosmic GRB. 



5. Discussion 



We have computed the mass loss rate in a magnetically driven wind emitted by a 
disk orbiting a stellar mass black hole. Detailed results are given for the case where 
the disk is optically thick to neutrinos but we have also obtained an approximate 
analytical expression, valid for any heating mechanism. From the mass loss rate the 
terminal Lorentz factor of the wind can be estimated. Large values of F require that 
severe constraints on the field geometry and the amount of heating below the sonic 
point should be satisfied. Another potential problem which was not addressed in the 
present paper is that relativistic MHD winds, at least in their simplest version fl3|| , 
can be quite inefficient in transferring magnetic into kinetic energy. 

An optimistic view of the situation would be to consider that the difficulty to 
obtain high Lorentz factors could just be a way to explain the apparent discrepancy 
between the birthrate of the sources in the hypernova scenario ~ 10~ 3 yr _1 /galaxy 
(for very massive stars) and the observed GRB rate < 10~ 6 yr _1 /galaxy. Beaming 
alone cannot account for the difference which implies that the collapse of a massive 
star most generally fails to give a GRB. 

The pessimistic view naturally consists to conclude that the baryonic load of the 
wind emitted by the disk is so large that the Lorentz factor can never reach values of 
10 2 and more. The next step is then to rely on the Blandford-Znajek mechanism to 
produce the relativistic wind One should be careful however that magnetic field 
lines coming from the disk and trapped by the black hole will also carry frozen in 
wind material leading to a "contamination" of the Blandford-Znajek process by the 
disk. A possible loophole could be that the accretion time scale to the black hole is so 
short that the stationary wind solutions we have obtained are not valid. The amount 
of material extracted from the disk could therefore be much smaller still allowing a 
highly relativistic outflow to develop. 
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